New perspectives on the Ising model 
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The Ising model, in presence of an external magnetic field, is isomorphic to a model of localized in- 
teracting particles satisfying the Fermi statistics. By using this isomorphism, we construct a general 
solution of the Ising model which holds for any dimensionality of the system. The Hamiltonian of 
the model is solved in terms of a complete finite set of eigenoperators and eigenvalues. The Green's 
function and the correlation functions of the fermionic model are exactly known and are expressed 
in terms of a finite small number of parameters that have to be self-consistently determined. By 
using the equation of the motion method, we derive a set of equations which connect different spin 
correlation functions. The scheme that emerges is that it is possible to describe the Ising model from 
a unified point of view where all the properties are connected to a small number of local parameters, 
and where the critical behavior is controlled by the energy scales fixed by the eigenvalues of the 
Hamiltonian. By using algebra and symmetry considerations, we calculate the self-consistent pa- 
rameters for the one-dimensional case. All the properties of the system are calculated and obviously 
agree with the exact results reported in the literature. 



I. INTRODUCTION 



It is really very hard to say something new on the 
Ising model. The model, originally proposed by Lenai in 
1920, was exactly solved for the case of an infinite chain 
by Ising£ in 1925. Since then, thousand and thousand 
of articles and several books have been published on the 
subject. The reason is that the model is very simple, but 
still can be considered as the prototype for systems sub- 
ject to second order phase transitions and can be effec- 
tively used for studying critical phenomena. Moreover, 
the model and its generalizations and modifications can 
also be used for studying a large variety of physical sys- 
tems. We do not attempt to summarize the enormous 
work done in these 80 years; it would go well beyond the 
purpose of this article. An excellent historical presenta- 
tion of the Ising model can be found in Ref.0, although it 
is old and obviously not updated. With no pretension of 
being exhaustive and complete, we here summarize the 
principal approaches used in these 80 years. 

A basic tool is the transfer matrix method-' 5,6,7 . By 
means of this approach, Onsager 8 in 1944 succeed to 
give an exact solution of the model for a cubic two- 
dimensional lattice in absence of external magnetic field. 
The theory of spinors and Lie algebra was used to sim- 
plify the Onsager solution&iS. Among the exact re- 
sults for the two-dimensional case, the calculation of 
the magnetization^! and the writing of the spin corre- 
lation function in the form of a Toeplitz determinant 12 
have to be mentioned. Other simplifications of the 
Onsager solution have been obtained by means of 
the Jordan- Wigner transformation and fermionization 
methodsi^ii^iiii. Different approaches are based on com- 
binatorial methoda 16 i 17 i 18 and pfaffian methoda 19 i 20 i 21 i 22 . 
More recent approaches have seen the Ising Hamiltonian 
expressed as a Gaussian Grasmannian action 2 ^ 2 *. Along 
this line, use of operatorial symmetries that simplify the 
algebra of the transfer matrix has led to the calculation 



of the partition function for a large class of latticesS 5 ^ 2 ^. 

Many approximation methods have been used with the 
goal of obtaining an expression for the partition function 
valid over a large temperature range: mean field theory, 
Bethe approximation 27 , cluster variational methods^, 
Monte Carlo simulations, series expansions. The spin 
correlation functions have been studied at the critical 
temperature^ 9 , and in the asymptotic regionSS* 3 ^. To 
study critical phenomena and critical indices, tools like 
series expansions 3 ^ 3 - 3 ^, scaling^SLS, renormalization 
group theor y 38 ! 39 have been used. 

In spite of the tremendous work done, many problems 
remain unsolved. The exact partition function in a fi- 
nite magnetic field is still unknown for dimensions higher 
than one. Very few exact results have been obtained for 
the three-dimensional model. There is no exact solution 
for the two-layer Ising model either. Most of all, a gen- 
eral approach which works in all dimensions and under 
general boundary conditions, although in some approxi- 
mation, is needed. 

In a recent work—, we have shown that there is a 
large class of models which are exactly solvable in terms 
of a finite number of parameters that have to be self- 
consistently calculated. The purpose of the present paper 
is to apply the method proposed in Ref. |4(J to the Ising 
model and to show that an exact solution of the model 
does exist for any dimension. In Section 2, we introduce 
the Ising model for a d -dimensional cubic lattice and 
show that the model is isomorphic to a system of localized 
spinless interacting particles, satisfying the Fermi statis- 
tics. In Section 3, the Hamiltonian of the latter model is 
solved, that is, a complete finite set of eigenoperators and 
the relative eigenvalues are determined. Then, as shown 
in Section 4, the exact form of the retarded Green's func- 
tion (GF) and of the correlation function (CF) can be 
obtained. In Section 5, we derive a set of equations for 
determining the charge/spin correlation functions. As 
the composite operators do not satisfy a canonical al- 
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gebra, the GF, the CF and the charge/spin correlation 
functions depend on a set of internal parameters not cal- 
culable by the dynamics. For the one-dimensional case, 
by means of the composite operator method^i*22i2ii, we 
calculate these internal parameters (Section 6) and the 
charge/spin correlation functions (Section 7). Although 
obvious, it is worth noticing that all the results reproduce 
the exact solution known in the literature. 

What are the advantages of the present method and 
what is new in the context of the Ising model? We present 
a new scheme of calculations for treating the model. The 
scheme is general and can be applied to any dimension. 
In the framework of this scheme we show that the model 
is always solvable for all dimensions. The energy spectra 
of the system are known. In the one dimensional case we 
show that the energy scales determined by the spectra 
rule the behavior at the critical temperature. It is rea- 
sonable to expect that this is true also for higher dimen- 
sions. General relations among different spin correlation 
functions have been obtained. These are exact relations 
and might be used to check the consistency of some ap- 
proximate treatments or numerical calculations. In order 
to get quantitative results for the cases of two and three 
dimensions we have to determine a finite small number 
of parameters. All the properties of the Ising model, the 
magnetization, the thermodynamical quantities, the spin 
correlation functions, depend on these parameters that 
have to be self-consistently determined. By using algebra 
and symmetry considerations we calculate these param- 
eters for the case d = 1. Extension of the calculations to 
higher dimensions is under investigation. 

II. THE ISING MODEL 

The Ising model, in presence of an uniform external 
magnetic field h, is described by the following Hamilto- 
nian 

H Ismg = £ W)S(i) - hJ2S(i) (1) 

ij i 

S(i) are spin variables, residing on a <i-dimensional Bra- 
vais lattice of N sites spanned by the vectors Rj = i . 
The variables S(i) takes only two values: up or down, 
or more simply S(i) = ±1. For a hypercubic lattice of 
lattice constant a with nearest neighbor interactions, the 
exchange matrix Jy is given by 

Jy = -2dJay ay = i E e ik (R *- R ^a(k) 

d k (2) 

a ( k ) = 3 E cos(k„a) 

n=l 

where d is the dimensionality of the system and k runs 
over the vectors in the first Brillouin zone. The exchange 
constant J can be positive or negative, and accordingly 
the coupling will be ferromagnetic or antiferromagnetic. 
According to J2J, the Hamiltonian Q can be rewritten 



as 

H Ismg = -dJ^ S(i)S a (i) - hJ2S(i) (3) 

i i 

where 

s«(i) = X>aS(i) (4) 
j 

It is worth to recall that the Ising Hamiltonian is 
invariant under the transformation 

h-*-h (5) 

Let us consider a system of N e interacting spinlcss 
fermions residing on the same lattice and let c(i) and 
(y{i) be the related annihilation and creation operators. 
These operators are Heisenberg fields [i — (i, i)] satisfy- 
ing canonical anticommutation relations 

{c(i,t),ct(j,t)} = 5 y 

{c(i, t), c(j, *)} = {ct (i, £), ct (j, t)} = W 

As a consequence of the algebra ©, each site can be 
occupied at most by a single particle. The occupation 
number of the site i, v{i) = c'(i)c(i), takes only the values 
and 1. By taking into account two-body interactions, 
the Hamiltonian for such a system reads as 

H =-E mo + \ E nuMOKi) (7) 

i ij 

where /i is the chemical potential and V(i,j) is the po- 
tential. This model Hamiltonian can be connected to the 
Ising model by defining 

v(i) = ±[l + S(i)} (8) 

It is clear that 

!/(*) = 5(i) = -l rq , 

v{i) = 1 & S(i) = +1 w 

By substituting (JSJ into (J7J) and by considering only a 
nearest-neighbor potential we can rewrite the Hamilto- 
nian Q in the following form 

H = E -hY,S{i)-djY,S(i)S a (i) (10) 

i i 

where we defined 

E = (-y+\Vd)N 

h=\(n-Vd) (11) 

J=-\Vd 

Hamiltonian (|10fl is just the Ising Hamiltonian (J3J) as we 
have the equivalence 

Hlsing = H — Eq (12) 
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The relation between the partition functions is 

Z H = e- 0Ea Z Ismg (13) 

Then, the thermal average of any operator A assumes the 
same value in both models 



{A{y)) H = (A(S)) Is 



(14) 



According to this, we can choose to study either one or 
the other model and get both solutions at once. We de- 
cide to put attention to the model Hamiltonian (JJJ , which 
for a nearest-neighbor potential reads as 

H = -n^v{C) + Vd^v{i)v a {C) (15) 



that for p > 1 the number operator v(i) = c*(i)c(i) sat- 
isfies the following algebra 



u(i)]P = [J{i)c(i)] p = v(i) 



(22) 



Therefore, the hierarchy of composite operators i|2U|) 
must close for a certain value of p and we should be 
able to derive a finite closed set of eigenoperators of the 
Hamiltonian. To this purpose, on the basis of (|22(l the 
following fundamental property of the field [^ a (z)] p can 
be established 



2c/ 



(23) 



where 



(16) 



The spin-inversion symmetry |J5j of the Ising model 
© corresponds to the particle-hole symmetry exhibited 
by the Hamiltonian (|15fl . In particular, we have that the 
chemical potential as a function of v — (v(i)) scales as 



III. 



fi(l -v) = 2dV - (i(v) 



COMPOSITE OPERATORS AND 
EQUATIONS OF MOTION 



(17) 



It is immediate to see that the charge density operator 
v(i) satisfies the equation of motion 



(18) 



Then, standard methods based on the use of equations of 
motion and Green's function (GF) formalism are not im- 
mediately applicable. Indeed, it is easy to check that the 
causal propagator (T[v(i)v(j)]) [T is the chronological 
operator] and the correlation function (v(i)v(j)) assume 
the form 

(T[u(i)v(j)]) = (K0"0')> = ^^> ik -( R *- R ^r(k) 

where T(k) is the zero frequency function^ which cannot 
be calculated by means of the dynamics^. 

Then, in order to solve the Hamiltonian Ijl5|) let us 
consider the composite operator 

^(i) = c(0KWf- 1 {p=M } (20) 

This field satisfies the equation of motion 

iJ^pW = bP P (i),H] = -^j p (i) + 2dV^ p+1 (i) (21) 

By taking higher-order time derivatives we generate a 
hierarchy of composite operators. However, we observe 



where the coefficients A^m satisfy the relation 

2<7 



(24) 



m— 1 



The proof of Eq. I12MI) and the explicit expressions of the 

coefficients Am are given in Appendix A for the cases 
d = 1,2,3. We now define the composite operator 



Mi) 



c(i)v a (i) 



(25) 



After (J2HJ), this field is an eigenoperator of the Hamilto- 
nian (flT|) 

i^ (d \i) = [^ {d) ii) 1 H]=e {d) M ) {i) (26) 

where the (2d + 1) x (2d + 1) matrix e^ d \ the energy 
matrix, is defined in Appendix B. It is easy to see that 
the eigenvalues En of the energy matrix are given by 



= -n + (n-l)V 



l,2,..-,2d+l (27) 



The Hamiltonian l|15fl has been solved since we know 
a complete set of eigenoperators and eigenvalues, and 
we can proceed to the calculations of observable quanti- 
ties. This will be done in the next Sections by using the 
Green's function formalism . 



IV. RETARDED AND CORRELATION 
FUNCTIONS 

We define now the thermal retarded Green's function 
G&(i,j) = (R[^(i)^(j)}] 

= e^-t^U^^^Hj)}) (28) 
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where (■ • ■ ) denotes the quantum-statistical average over 
the grand canonical ensemble. By introducing the 
Fourier transform 

1 



N ^ (2?r) 



J dwe ik - (Rl - Rj) - iw(t ^ tj) G (d) (k,w) (29) 



and by means of the Heisenberg equation (|26[) we obtain 
the equation 



[u-£W]GM(k,w)=j( J '(k) 



(30) 



where /^ d ^(k) is the Fourier transform of the normaliza- 
tion matrix, defined as 



(31) 



= ({^(i,t),^)t(j,i)} 

= J_^ e ik-(R l -R 3 ) / (<i)( k ) 



The solution of Eq. (gO} is 



2d+l 



(d.n) 



(k) 



n w-£;.i d) +i(5 



(32) 



n=l 



The spectral density matrices cr^i (k) are calculated by 
means of the formula— » 



«r2- w J(k) = nw^[nw ] -x i g) (k) 



(33) 



where fi^ d 'is the (2d 4- 1) x (2d+ 1) matrix whose columns 
are the eigenvectors of the matrix e^ d \ The explicit ex- 
pressions of O^are given in Appendix B. The spectral 
density matrices cr( d: ")(k) satisfy the sum rule 



2d+l 



53 [Ei d) ] p a^' n \k) = M^\k) 



(34) 



where M( d ' p )(k) are the spectral moments defined as 
M Wp)( k ) = F T /{{idjdtf ip (d \i, t), ip {d) Hj, f)} 



(35) 

F.T. stays for the Fourier transform. It is a consequence 
of the theorem proved in Ref. |44] [see also pag. 572 in 
Ref. l4l| that the spectral density matrices, for d = 1, 2, 3, 
satisfy the sum rule pi|l. The explicit expressions of 
J (d) (k) and cr (d ' n) (k) are given in Appendices C and D, 
respectively, for the cases d — 1,2,3. The correlation 
function 



(36) 



can be immediately calculated from (|32|l by using the 
spectral theorem and one obtains 

+ 00 

= lyJL f dwe i(R,-R J )-iu»(t 1 -t,) C ,(«o/ k ^ j (37) 

N ^ (2tt) 7 v > / v ; 



2d+l 



<7 (d) (k,w) = 7T ^ S[lj - £jW]TW<7^' n) (k) (38) 



with 



T{ d) = 1 + tanh 



E. 



(d) 



2k B T 



(39) 



Equations l|32|) and l|38|l are an exact solution of the 
model Hamiltonian (JTSJ. One is able to obtain an 
exact solution as the composite operators ip p (i) — 
c(i)[v a (i)] p ~ 1 constitute a closed set of eigenoperators 
of the Hamiltonian. However, as stressed in Ref. EH the 
knowledge of the GF is not fully achieved yet. The alge- 
bra of the field ip( d \i) is not canonical: as a consequence, 
the normalization matrix 1^ (k) in the equation ( I3U|) 
contains some unknown static correlation functions, cor- 
relators (see Appendix C for explicit calculations), that 
have to be self-consistently calculated. According to the 
scheme of calculations proposed by the composite oper- 
ator method^i^24^(COM), one way of calculating these 
unknown correlators is by specifying the representation 
where the GF are realized. The knowledge of the Hamil- 
tonian and of the operatorial algebra is not sufficient to 
completely determine the GF. The GF refer to a spe- 
cific representation (i.e., to a specific choice of the Hilbcrt 
space) and this information must be supplied to the equa- 
tions of motion that alone are not sufficient to completely 
determine the GF. Usually, the use of composite oper- 
ators leads to an enlargement of the Hilbert space by 
the inclusion of some unphysical states. Since the GF 
depend on the unknown correlators, it is clear that the 
value of these parameters and the representation are in- 
timately related. The procedure is the following. We 
set up some requirements on the representation and de- 
termine the correlators in order that these conditions be 
satisfied. From the algebra it is possible to derive several 
relations among the operators. We will call algebra con- 
straints (AC) all possible relations among the operators 
dictated by the algebra. This set of relations valid at 
microscopic level must be satisfied also at macroscopic 
level, when expectations values are considered. Use of 
these considerations leads to some self-consistent equa- 
tions which will be used to fix the unknown correlators 
appearing in the normalization matrix. An immediate 
set of rules is given by the equation 



<- u 'Hnc«»Hi)) =^E^ / <^"' ,:, (k^! (40! 
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where the l.h.s. is fixed by the AC and the boundary 
conditions compatible with the phase under investiga- 
tion, while in the r.h.s. the correlation function C(k, ui) 
is computed by means of equation of motion [cfr. Eq. 

Another important set of AC can be derived by observ- 
ing that there exist some operators, O, which project out 
of the Hamiltonian a reduced part 

OH = OH a (41) 

When Hq and Hi = H — Hq commute, the quantum 
statistical average of the operator O over the complete 
Hamiltonian H must coincide with the average over the 
reduced Hamiltonian Hq 

Tr{Oe- 0H } = Tr{Oe-^ H "} (42) 

Another relation is the requirement of time transla- 
tional invariancc which leads to the condition that the 
spectral moments, defined by Eq. 1)351 ) . must satisfy the 
following relation 

M^(k) = [M b ( a d ' p) (k)]* (43) 

It can be shown that if 1|43[) is violated, then states with 
a negative norm appear in the Hilbert space. Of course 
the above rules are not exhaustive and more conditions 
might be needed. 

According to the calculations given in appendices 
C and D, the GF and the correlation functions de- 
pend on the following parameters: external parameters 
(fj,,T,V), internal parameters (C^™, C^j™, •• • C^ld )> 
and (k*- 1 -*, ■ ■ ■ k^ 2 ^), defined as 

ClT = (< )Q (*)^ d)t «) (44) 



«W = ([v a (i)] p ) (45) 

The parameters Cjfy a are determined by means of their 
own definitions (|44|) . where the r.h.s. is calculated by 
means of ( I37|) - (|38|l . This equation gives 

, 2d+l 

C ( d)a = £ T W-5>(k)^<«)(k) (46) 

n=l k 

From the results given in the Appendices C and D, we 
see that the spectral density matrices have the form 

CT K«)(k)=A^ n) +a(k)A^' n) (47) 

where the matrices A and Ai do not depend on momen- 
tum k. Putting 147|) into i|46|) we obtain 

1 2d+l 

° {d)a = J7i £ 4 d) Ai d,n) (48) 

n— 1 



Calculations given in the Appendices C and D show that 
the matrices A^f' 71 ^ are linear combinations of the matrix 
elements cj d p )a . Then, Eq. gHJ gives a system of homo- 
geneous linear equations. The determinant of this sys- 
tem is only function of the external parameters /i, T, V. 
This function will vanish only if there is a particular re- 
lation among these parameters. Since these parameters 
are independent variables the only solution is that all the 
matrix elements must vanish 

C[ d) p a = (49) 

The matrices A^'™^ are zero and the correlation func- 
tion C^(k, uo) does not depend on momentum, as we 
expected. In the coordinate space the CF takes the ex- 
pression 

, 2d+l 

CV\i,j)=8 i ±Y J T n A a d ^e-^ ) ^ (50) 

n=l 

The correlation function depends on 2d internal parame- 
ters: kS 1 ', • • • , n( 2d \ In order to determine these param- 
eters, we use the Pauli condition (|40|l which gives the 
self-consistent equations 

K W-AW=Cg +1 (p = 0,l,-.-2d) (51) 

where C[ d p +1 = \^[ d \i)'>pp d ^ is calculated by means 
of (|5t)(l . New correlation functions 

A<*> - {v{i)[v a {i)] p ) (52) 

appear and the set of self-consistent equations (|51|l is 
not sufficient to determine all unknown parameters. One 
needs more conditions. In the case of one-dimensional 
systems these extra conditions can be obtained by using 
the property (g3J). 

V. CHARGE CORRELATIONS FUNCTIONS 

In Sections 3 and 4, we have solved the problem of the 
Ising model in terms of a set of local parameters, defined 
by l|45l) and (|52ll . In this Section, we want to show how 
we can calculate non-local correlation functions. Let us 
define the causal Green's function (for simplicity in this 
Section we drop the superindex (d)) 

F c {i,l,j) = (TW)tf(l)]vV)) 

-flfe-tiK^OV^Mj)) (53) 

the retarded and advanced functions 

F R > A (i,l,j) = (R,AW)tf{l)]vV)) 

= ±6[±(t i -t l )]({i>(i),^(l)}v(j))(5A) 
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the correlation functions 



DW(i,l,j) = (ip(i)^(l)v(j)) 
D^(i,l,j) = (^(l)i>(i)v(j)) 



(55) 



where ip(i) is the composite field defined in (|25l) and we 
used the fact the field operator v{j) does not depend on 
time. The Fourier transforms of these quantities read as 



(i, J, j) = ±J diue-^-^D^ (i, 1, j; w) 1 j 

where Q = C,R, A. By means of the equation of motion 
(l26l) we have 



(w-e)i*(i,l,j;w) = J(i,lJ) 



(w - 6)2^(1, 1, j;w) = 



(57) 



(58) 



where the matrix J(i, 1, j) is defined as 

J(i,l,j) = <Wi,t),^t(i, t )}„(j)) (59) 

The most general solution of Eq. (|57|) is 

~ U -E n 

n—1 

-m5{uo - E n )g^ n \iX])] (60) 

where 

2d+l 

r^(i, 1, j) - tt a „ n-i J c6 (i, 1, j) (61) 

c=l 

while the function gQ( n '> (i, 1, j) must be determined. P 
denotes the principal value. By recalling the retarded 
and advanced nature of F R ' A (i,l,j), it is immediate to 
see that 

g R ^(i,l,j) = -g A ^(i,l,j) = rW(i,l,j) (62) 
Therefore 

n—1 

The solution of JSSJl is 



2d+l 



D^ t (i,l,j;w)= E %-^n)d** t(n) (i.lJ) 



71=1 

2d+l 



D^(i,l,j;w) = E %-^)i* t,lw (i,lj) 



(64) 



where the matrices ef^^i, I, j) and d^^ n \\, 1, j) have 
to be determined. From the definitions l|53|) - (|55|l we can 
derive the following exact relations 

F R (i,l,j)+F A (i,l,j) = 2F c (i,l,j) - (bP(i),^(l)Hj)) 
F R (i,l,j) - F A (i,l,j) = ({^(0,^(0 Mi)) 

(65) 

A relation between the two correlation functions 
DW\i,l,j) and D^^(i,l,j) can be established by 
means of trace properties. Indeed, it is straightforward 
to derive a KMS-like relation 

(^(l)il>(i)v(j)) = (i>(i,ti-i{3)^(l)is(j)) 

+5 li (iP(i,t i -ip)^(l)) (66) 

By recalling the definitions l|55[) . this last equations can 
be written as 

D**+(i, 1, j; w) = e~^[D^ (i, 1, j; w) + 1; w)] 

(67) 

where C(i, l;w) is the fermionic correlation function 
[see Eq. ( I37[) - H38fl ]. Therefore, the anticommutator 
\{ip(i),ip'(J)}v(J)) in l|65(l can be expressed in terms of 
the correlation functions as 

x[(l + e-^)Z>^ f (i,l,j;o;) + 5y e -^C(i,l;a;)] 1 J 

Analogous expression holds for the commutator. By 
means of and by recalling that [see Eqs. lpT7|) - (|5"8"|) ] 



(w-e)C(i,l;w) =0 

C(i,l;a;) = 2 E 1 5(a;-^ n ) C («)(i,l) ,<,<M 



we find that equations (|65|l have the following form 



2d+l 



E ^-P„){.9 C ^(MJ)- ^[(l-e-^) 

71=1 l' u 7 

xd** t W(i,IJ)-<y u e-' Sw c( n )(i,l)]} = 



2d+l 



E ^-^{rWaij)- ^[(l + e-^) 

71=1 



(71) 



By recalling that = —fx + (n — \)V 1 the solution of 
and (El) is: 

2tt 



d**t(»)(i,i,j) = __^_ T (»)(i,l,j) 

- C (")(i,l) (72) 



1 + e^" 



1 + e -^» 



3 C(n) (i,lJ) 



1 + e-^" 1 ' ' Jj 
2tt l + e-^» { , } (76) 
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By putting J7U and (JJSJ into (fSUJ) and H2} we have 

2d+l iB„(t ( -t,) 

^ (Ui) = ETT^r [rW(MJ) 

-^<5 U e-^«cW(i,l)] (74) 



2d +! e - e -/3-E„ 



+ e 



[r (n) (MJ) 



+i-,5 ljC (»)(i,l)] 



£ - )(UJ) 
^ i 

n=l 

1 P -/3£„ 



(75) 



Lij — E n +iS uj — E n —\5 

2 ^c(")(i,l)e-^ 



1 



1 



cj — i? n +i<5 cj — E n —i8 
From the study of the fermionic sector we have 

2tt 



,(' 



) (M) = 



1 + e-( iE - 



-<W (ri) 



(76) 



(77) 



where er^are the spectral functions given in Appendix 
D. Then, at equal time l|74j) becomes 



2d+l 



,-PE n 



x[rW(i,IJ)-M U 1 + ^ g „ ^ (n) ] (78) 

The system l(78j) gives a system of linear equations 
for the quantities Z?^ (i, l,j). Since the inhomogeneous 
terms in this system are proportional to Sn, it is clear 
that (i,I,j) oc 5 a . Then, we will take i = 1 and we 
write 

D(iJ) = (^(i)^(i) n (j)) (79) 
From l)78|l we have the system of equations 

2d+l 1 

^frW(i,j)-fe- 

+ e'- 



2d+l 1 1 



where 



, , 2d+l 

7-w(i,j) = n on g n-ijc6(i,j) 
j(i,j) = ^(i)/(i)l(j)) 



(80) 



(81) 



From its own definition l|79|) and by using the recurrence 
relation f 123(1. the matrix £>(i, j) has the following struc- 
ture. 

(i) One dimension 



. D hl D 1>2 D h3 
£>^(i,j) - ( £>1,2 Dl,3 #2,3 
-Dl.3 ^2,3 -D3.3 



(82) 



D hp (i, j) - J^ 1 ) (i, j) - A^Hi, j) p = 1, 2, 3 



A>,3(i, J) = E A^ +1) £>i >m+1 (i, j) p = 2, 3 

m— 1 



(ii) Two dimensions 



£ (2) (U) = 



/ £»!,! Di >2 Dl,3 £>1,4 #1,5 \ 

#1.2 £>1,3 £>1,4 #1,5 #2,5 1 

I) .:; Dl,4 £>1,5 #2,5 /A; ", 

Al.,4 -Dl.5 ^2,5 -D3.5 

V #1,5 #2,5 ^3,5 #4,5 £>5,5 / 



(83) 



(84) 



£>i, p (i, j) = A"**- 1 ) (i, j) - A(P-D(i, j) p = 1, 2, • • • , 5 
A>,5(iJ) = E ^ +3) D 1>Tn+1 (i, j) p = 2, 3, • ■ • , 5 

m— 1 



(iii) Three dimensions 



£ (3) (U) = 



D12 D 13 ■ ■ ■ D lj5 D 2 ,7 



Di : e -D1.7 • • ■ D 5 ^ D 6 ,7 

VD17 132,7 ••■ A>,7 ^7,7/ 



(85) 



(86) 



j) = (i, j) - A(p-D(i, j) p = 1, 2, • • • , 7 

A>,r(iJ) = E ^ +5) £>i,m+i(i, j) P = 2, 3, • ■ • , 7 

m=l 

(87) 

With the definitions 



ifW(iJ) =< [^ Q (i)]Mj) > 
AW(iJ) =< 2,(i)[*/*(i)]Mj) > 



Then, we only need to calculate the matrix elements 
Dx lP (i, j) (p = 1, 2, • • • 2d + 1). The matrix J(i, j) can 
be obtained from the normalization matrix J(i,j) = 
({^>(i, t),tp'(j, t)}), calculated in Appendix C, by means 
of the following substitution 



,(p) 



(89) 



Then, the matrices r^(i,j) have the same expressions 
of the spectral matrices when the following substi- 
tution 



I, 



ab 



Jab(h3) 



(90) 
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is made. It can be seen that for j = i and j = i Q the sys- 
tem (I8UII is exactly equivalent to the system l|51[) . Then, 
it is enough to consider the case j ^ i, i a . In this case, 
the system (I8U|) becomes 



where 



, 2d+l 

£>(i,j) = -^T n r(")(i,j) 

Z n=l 



(91) 



with T n given by i|39|) ■ The system <|91|) gives a set of ex- 
act relations among the correlation functions. We might 
think to solve this system by induction method, since 
some of the first correlation functions can be expressed 
in terms of the basic parameters and A^. However, 
when we do this, we immediately see that the number of 
equations is not sufficient to determine all the correlation 
functions and we need more equations. Once again, this 
can be done for the one dimensional system, as we shall 
see in the next Sections. 



VI. SELF-CONSISTENT EQUATIONS FOR 
ONE-DIMENSIONAL SYSTEMS 



Until now the analysis has been carried on in complete 
generality for a cubic lattice of d dimensions. We now 
consider one-dimensional systems, and in particular we 
will study an infinite chain in the homogeneous phase. 
For simplicity of notation we shall drop the superindex 
(d). By means of HC6|) and l)Dl|) - (|D2j) the set of equations 
B51J1 gives the linear system 

Ti-2 + (2-3Ti+4T 2 -T 3 )^ 
+2(T 1 -2T 2 + r 3 ) K ( 2 ) =0 

(2T 2 - T 3 - 2)v - 2(T 2 - T 3 )^ 2 ) + 2A« = 1 > 
(T 2 - T 3 )v - (2 +T 2 - 2T 3 )k^ + 2A< 2 ) = 

where, because of translational invariance, we put 

v = = (93) 

It is immediate to see that for fi = V, the solution of the 
first equation in 192JI for T > is 



1 

v = — 
2 



for ji = V 



(94) 



This is in agreement with the particle-hole symmetry en- 
joyed by the model [see l(T7|> ]. Recalling JSJ) and (jTT|) . 

this situation corresponds to the zero magnetization of 
the Ising model in absence of external magnetic field. 
Coming back to general value of /i, it is clear that Eqs. 
(|92|l are not sufficient to specify completely the 4 param- 
eters v, k( 2 \ A^ 1 - 1 , A^ 2 ) and we need another equation. A 
fourth equation can be easily obtained by means of the 
algebra. We observe that 



c\i)v(i) = 

This relation leads to 



(95) 



(96) 



H = H -2Vv(i)v a {i) 



(97) 



By means of the requirement (|42|l the correlation function 
C\.k — (c{i)v {i)[v a (i)\ k \ can be written as 



where 



(98) 



(99) 



and (• • • ) denotes the thermal average with respect to 
Hq. Let us define the retarded GF 

G<5(i-0 = (R[c(i,t)J(i,t')][v a (i)] k -\ 

= JL f dwe-M^Gg)^) (ioo) 



By means of the equation of motion 

[c(i),H ] = -fic(i) 

we have 



(101) 



(102) 



Recalling the relation between retarded and correlation 
functions, from l|102|) we obtain 



(103) 



By putting this result into l|98|l we have 

C hk = Ci,i <K(i)]^ X > (104) 
By noting that [^ Q («)] 2 can be expressed as [cfr. (|A4h ] 

[v a (i)? = lW a (i) + v(h>(i2)\ (105) 
we obtain from l|104|l the relations 

Ci, 2 = &,! ([u a (i)}) (106) 

Ci,3 = \[Cx,2 + C 1A {v{h)v{i 2 )) ] (107) 

Now, we observe^ that Hq describes a system where 
the original lattice is divided in two disconnected sub- 
lattices (the chains to the left and to the right of the site 
i). Then, in Hq -representation, the correlation function 
which relates sites belonging to different sublattices can 
be decoupled: 



(a(j)b(m)) = (a(j)) (6(m)) 



(108) 
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for j and m belonging to different sublattices. By using 
this property, invariance of Hq under axis reflection and 
(jl ()6|l we can write 



(v(ii)v(i 2 )) = {v{ii)) Q {v(i 2 )) 



Ci, 



(109) 



By putting (|109|> into (|107|l . we obtain the following self- 
consistent equation among the correlation functions 



^1,3 — oCl,2 



Cl,2 



(110) 



By means of (|37[1 and 1|38|) and the results given in Ap- 
pendices C and D, Eq. (|110fl takes the expression 



(4T 2 - 'iT^v 2 + [TiT 3 - 8 K ( 2 )(T| - TiT 3 )]i/ 
+2k( 2 )[2k( 2 )(T| - TiT 3 ) - T1T3] = 



(111) 



This equation together with equations (|9^|) gives a sys- 
tem of 4 self-consistent equations for the 4 parameters 
v, , A^ 1 ) , A^ 2 ) as functions of fi, T, V. By solving the 
set of linear equations l|9"2l with respect to , A^ , A^ 2 ^ 
as functions of v, we have 



.(2) 



A< 2 > = 



2 - Ti + i/(-2 + 37\ - 4T 2 + T 3 



2(Ti 
1 



2T 2 + T 3 ) 

-{ (2-r 1 )(r 2 -r 3 ) 



(112) 
(113) 



2(Ti-2T 2 + r 3 ) 
+i/[Ti(2 + T 2 - 2T 3 ) + T 2 (T 3 - 6) + 4T 3 ]} 

4(7V^T^) {(2 - Tl)(2 + T2 - 2T3) (U4) 
+u[-4 - 1QT 2 + T 1 {6 + T 2 - 4T 3 ) + 6T 3 + 3T 2 T 3 ]} 



To calculate the parameter v let us put i|112[l into l|lll|) 
and solve with respect to v. We have two roots. One 
solution corresponds to an unstable state with negative 
compressibility and must be disregarded. By picking up 
the right root, and by using the relation 



T 3 = 



2T 2 2 (2-T 1 ) 



T 1 (2-T 2 ) 2 + T 2 (2-T 1 ) 



we find 



1 

V= 2 



1 + 




Ti 



Ti - 2TiT 2 + 2T 2 



(115) 



(116) 



As shown in Appendix E, the solutions (|112fl - i|116[l ex- 
actly correspond to the well-known solution of the ID 
Ising model, obtained by means of the transfer matrix 
method. We could manipulate the expression (|116Jl and 
the ones for iv- 2 \ A^ and , obtained by substituting 
(|116fl into l|112|) - (|114fl . m order to reproduce the expres- 
sions of the Ising model, given in Appendix E. However, 
we prefer to maintain the present expressions as the fol- 
lowing discussion will be more transparent. In Section 
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FIG. 1: The particle density v is plotted as a function of the 
chemical potential at various temperatures for V = — 1 (top) 
and V = 1 (bottom). 



3, we have seen that in the present model, in the one- 
dimensional case, there are three energy scales 



Ei = -11 

E 2 = -n + V 

E 3 = -(i+ 2V 



(117) 



At zero temperature, the three functions Ti, T 2 and 
T 3 are not analytical functions at the points /1 = 0, = 
V, and fi = 21^, respectively, and we expect that the 
parameters i/, k 

(2) )A (i) )A (2) exhibit 

some discontinuous 
behavior at these points. As shown in Fig.^ in the limit 
T — > the particle density v has a discontinuity at ji = V 
for the case of negative V (i.e. J > 0, ferromagnetic 
coupling) and two discontinuities at (i = and /i = 2V 
for the case of positive V (i.e. J < 0, antiferromagnetic 
coupling). Here and in the following, we take \V\ = 1: 
all energies are measured in units of \V\. In particular, 
the particle density increases by increasing fj, from zero 
to one. At zero temperature, in the ferromagnetic case 
v is zero for fi < — \V\ and equal to one for fi > — |V|; 
in the antiferromagnetic case v is zero for /i < 0, jumps 
to 1/2 and exhibits a plateau, centered at fi = V, in the 
region < \i < 2V, jumps to the value 1 for \l > 2V. 
The parameter has a behavior similar to v. 



10 



1.2 p 

1.0 - 

0.8 - 

0.6 - 

0.4 - 

0.2 - 

0.0 - 
-0.2 



1.2 
1.0 
0.8 
0.6 
0.4 
0.2 
0.0 
-0.2 



V=-l 



/ 

/ 

f ' 
§ 



T=0 
T=0.25 
T=0.5 
T=1.0 



6 


-4 -2 





2 4 ( 






li 














T=0 








T=0.1 




• _ - - ' 




T=0.5 








-T=1.0 












1 / 




V=l 




/'/' 








/ 

/ _ 






/[ 






























M 



FIG. 2: The correlation function A' 1 ^ is plotted as a function 
of the chemical potential at various temperatures for V = — 1 
(top) and V = 1 (bottom). 



>( K ( 2 )-A( 2 ))^(^-AW)+ ( ' y ; A(1))2 (1 



1-1/ 



19) 



The Ising model in one dimension can be described in 
terms of only two parameters: v and A^ 1 ). 



VII. CHARGE CORRELATION FUNCTIONS 
FOR ONE-DIMENSIONAL SYSTEMS 

The system of equations (|91ll establishes some 
relations among the non-local charge correla- 
tion functions K^(i,j) =< [v a (i)] p v(j) > and 
A (p) (iJ) =< v(i)[v a (i)] p v(j) >. As already discussed, 
the number of equations is not sufficient to determine 
completely the charge correlation functions, and one 
needs more equations to close the system. In the 
one-dimensional case a fourth equation can be eas- 
ily obtained by algebraic considerations. Recalling 
that ^(ije^ 1311 = c^(i)e~P H ° we can easily derive the 
following result 

A«(i, j) = K^^j) - C M (K«]Mj)) (120) 

Now, for j > i + a (because of invariance under axis 
reflection we could choose j < i — a as well), using the 
property l|108|) 



x(Kj))o+2^ + «)^(j))o( 121 ) 



In Fig. |21 we give the parameter A^ as a function of 
For the ferromagnetic case the behavior is similar to 
that of v . Instead, in the antiferromagnetic case at 
T = 0, exhibits only one discontinuity point at /i = 2V, 
where jumps from zero to one. The parameter A' 2 ' has 
a behavior similar to A^. The different behavior ex- 
hibited by the pairs (v, k (2) ) and (A (1) , A^ 2 )) for V < is 
naturally due to the antiferromagnetic correlations, when 
we recall that the pair (A^, A^ 2 )) describes a correlation 
between two first neighboring sites, while describes 
correlations between two second neighboring sites. Of 
course in the point of discontinuity the two limits T — > 
and [L — > fi c are not interchangeable. As we shall see in 
the next Section, the 4 local parameters v, kP\ A^ 2 -* 
are really basic since all the properties of the model are 
described in terms of them. It is worthwhile to note that 
some simple relations can be established among the pa- 
rameters 



+ a)u{i - a)v{j)) Q = {v a {i)) + a)u(j)) 

(122) 

Therefore, from 1120|) 

A^(i,j)=K^(i,j)-C 1 , 1 (u(j)) (123) 



A«(z,j) = K^{i,j)-ho x , 2 {v{j)) Q 

-±C hl (v(i + a)v(j)) (124) 



A< 2 )(m) = KW(i,j)-±[KU(i,j)-AW(i,j)\ 



(125) 



where we used Eq. IjlObj) ■ Equations (|123|) and l|124|l give 



2 K W 



(118) 



K?))o = 7^ ( °Hm)-A (0) (m)] (126) 
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C*12 



[K^(i,j)-M°\i,j)} (127) 



Putting i|127|l into (|125|1 we get the fourth self-consistent 
equation 

A^(t,j) = K^^j) 

-[\ + ^][K {1) {U)-^\i,3)] 

+^[^ (0) (i,i)-A^(i I i)] (128) 

It is straightforward to verify that (|128fl is identically sat- 
isfied for j = i, while for j = i a coincides with equation 
(TUT)]) . 

By adding (|128f) to l|91l) , we have the following system 
of equations for the non-local spin correlation functions 
K&(i,j) and AW(i,j) 

2A^ (i,j)-a K^(i,j)+ ai K^(i,j) = a 5 u (129) 



2AW (i, j) + a a #«(i,i) - 2a 3 ^ (2) (z,j) = (130) 



2A« (ij) + a 3 K^ (*, j) - ai K& (i,j) = (131) 



a( 2 )(m) + P o[# (1) (U)-a (1 Hm)] 

+p 1 A(°)(i,j)-tf( 2 )(i,j)=p 1 */ 
where we put 



(132) 



a = (37\ - 4T 2 + T 3 ) a 3 = (T 2 - T 3 ) 

ax = 2(Tx - 2T 2 + T 3 ) a 4 = (2 + T 2 - 2T 3 ) 

a 2 = (2T 2 - T 3 - 2) a 5 = (2 - T x ) 



Po 



Pi 



/l I Cl,2 ■ 

^2 Ci,i - 

°1,2 



(133) 



(134) 



and we used K^°'(i,j) = = v. To calculate for 

general |« — j'| we proceed by induction. 

Let us put i — j = ma and concentrate the attention 
on the spin CF A^°\m) = {v(m)u(0)). We start by ob- 
serving that 



A<°>(0) = i/ A(°)(1) = A( 1 ) 



(135) 



where the two parameters A^ have been calculated in 
the fermionic sector. By taking m = 2 we can calculate 
from the system 1(132 )1 that 



A(°)(2) = 2 K ( 2 )-!/ 



(136) 



The results (|135|) - (|136|1 an d use °f the relation l|118(l show 
that A(°)(m) for m — 0, 1, 2 can be cast in the form 



A(°)(m) = ^ 2 + j/(l-j/)p 
where the parameter p is defined as 

2/t (2)_„_„2 



i/ - i/- A« - j/ 2 



P : 



AW - i/ 2 



(137) 



(138) 



By using the expressions of the basic parameters given in 
Section 5, it is possible to check that |p| < 1. Then, we 
can introduce the Fourier transform and reexpress l|137|) 
as 



7r I a 



A^(m) = v 2 +Av{l-v) — / dk 



where 



42) 



.4 



27r J "" 1 + Bcos(fca) 

— n / a 



B 



(139) 



(140) 



Then, the CF A^(m + 1) can be calculated as 

7r j a 

f ika(m+l) 

A""(m + 1) = v 2 + Av(l-v)— \ dk— — 

V ' 1 2vr i l + Scos(fca) 

— n j a 

= v 2 + v{\ - v)p m+1 (141) 

Therefore, Eq. I|137fl is valid for any m. Recalling the 
definition of A^ ^ (m) , we can rewrite ()137ll under the form 



< v(m)v({) > -v 2 



= P 



(142) 



Also, from (|139fl we see that the zero frequency function 
[cfr. (|19fl ] has the expression 

T(fc) = v 2 (2n/a)6(k) + ^ ~ ' ^ (143) 
1 + B cos(ka) 

By putting the obtained expression of A^°\m) in Eqs. 
(|129l 'l- (ll32ll . we can solve the system. The solution gives 

K {1) {m) = v 2 + ~i/(l - */)(p m_1 +P m+1 ) (144) 



^{p m - 1 +p m+1 )-—f 
2ai a\ 



(145) 



AW( m ) = A (1 V + z/(l -v)x 



2a a 3 - aia 2 , m-1 



4a i 



2a 3 

2 

ai 



(146) 
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FIG. 3: The parameter p is plotted as a function of the chem- 
ical potential at various temperatures for V = — 1 (top) and 
V = 1 (bottom). 



-(p n - 

4ai 



04 
— P 

ai 



(147) 



In Fig. we give p as a function of \i for V = — 1 and 
V = 1 at various temperatures. We see that for negative 
1/ (ferromagnetic case), p is positive and various between 
zero and 1. For positive V (antiferromagnetic case), p is 
negative and various between —1 and zero. In particular, 
for negative V, p tends to 1 at /i = V in the limit T — > 0. 
Instead, for positive V, p tends to — 1 at fi = V in the 
limit T — > 0. This is seen in Fig. 0| where p is plotted 
versus T ai n = V = — 1 and at /i = V = 1. 

Let us now discuss the correlation functions. A^(m) 
is plotted against m for /i = V = — 1 [Fig. El (top)] and 
for (j, = V = 1 [Fig. (bottom)] at various temperatures. 
We see that when zero temperature is approached a long- 
range order of ferromagnetic and antiferromagnetic type 
is established, respectively. Also, we can see from ( I141fl - 
(I147H that for m — > oo and T ^ (i.e. not at the critical 
temperature) the spin correlation functions assume the 
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FIG. 4: The parameter p is plotted as a function of the tem- 
perature at \i = V = — 1 and n — V — 1. 
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FIG. 5: A (0) (m) is plotted against m for fj, = V = — 1 (top) 
and fj, = V = 1 (bottom) and various temperature. 



crgodic value. At the critical temperature T — we have 
breakdown of the ergodicity. 
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VIII. CONCLUSIONS 



The Ising model in presence of an external magnetic 
field is isomorphic to a model of localized spinless in- 
teracting particles, satisfying the Fermi statistics. The 
latter model belongs to a class of models always solv- 
able, as shown in Ref. |4fJ. On this basis, we have con- 
structed a general solution of the Ising model which holds 
for any dimensionality of the system. The Hamiltonian 
of the model has been solved in terms of a complete fi- 
nite set of eigenoperators and eigenvalues. The Green's 
function and the correlation functions of the fermionic 
model are exactly known and are expressed in terms of 
a finite small number of parameters that have to be self- 
consistently determined. By using the equation of the 
motion method, we have derived a set of equations which 
connect different spin correlation functions. The scheme 
that emerges is that it is possible to describe the Ising 
model from a unified point of view where all the proper- 
ties are connected to a small number of local parameters, 
and where the critical behavior is controlled by the en- 
ergy scales fixed by the eigenvalues of the Hamiltonian. 
The latter considerations have been proved from a quan- 
titative point of view in the one-dimensional case, where 
the equations which determine the self-consistent param- 
eters and the spin correlation function have been solved. 
For d = 1 all the properties of the system have been 
calculated and obviously agree with the exact results re- 
ported in the literature. Extension of the calculations to 
higher dimensions is under investigation. 



After the paper was completed and submitted for pub- 
lication, the author learned that approaches to the spin- 
1/2 Ising model based on a fermionization of the model 
have been previously reported in Refs. \^ and 0. The 
author wishes to thank the referee and Prof. L. De Cesare 
for putting these papers to his attention. In particular, in 
Ref. 0, Tyablikov and Fedyanin showed that the chain 
of equations for the double-time GF closes and the num- 
ber of equations is determined only by the co-ordination 
number, independently by the dimensionality This con- 
clusion agrees with the results given in Sections 3 and 
4. In order to close the set of equations for the fermionic 
correlation functions in the one-dimensional case, the au- 
thors of Ref. |4fj| assumed ergodicity and solved the sys- 
tem, obtaining the exact solution of the ID Ising model 
for an infinite chain. It should be remarked that ergod- 
icity breaks down for finite systems and at the critical 
points. Kalashnikov and Fradkin in Ref. used the 
spectral density method^ to derive a system of equa- 
tions for the correlation functions; also in this case the 
approach is valid for any dimension. However, the num- 
ber of equations is less than the number of correlation 
functions. 



APPENDIX A: ALGEBRAIC RELATIONS 

As mentioned in Section 3, the number density opera- 
tor v(i) — c^(i)c(i) satisfies the algebra 

u p (i) = v(i) p > 1 (Al) 

From this algebra an important relation can be derived, 
for the operator 

2d 

^« = E a y^') = Yd^ v{lm) (A2) 

j m=l 

where i m are the first neighbors of the site i. We shall 
discuss separately the cases of different dimensions. 

1. One dimension 

We start from the equation 

K«F = ^ E ( m ) ^ir-Vfer (A3) 

m=0 ^ ' 

After subtracting the terms m = and m = p, we can 
use the algebraic relation (|A1|) to obtain 

= l[2^)+a P K*iM*2)] (A4) 

with 

« P = E(rO =2P - 2 (A5) 

m=l V 7 

From (|A4|) , by putting p — 2 we obtain 

v(ii)v(i 2 ) = 2[v a (i)] 2 -V a (i) (A6) 

By substituting i|A6(l into (|A4|) we have the recurrence 
rule 

2 

[z/w = E A ™V«r (a?) 

m— 1 

where 

4 P) = ^(2-a p ) = 2 2 ^-l (A8) 
A P) = i2a p = 2(l-2 1 ^) (A9) 
We note that the coefficients 4™ satisfy the relation 

2 

E A ™ = 1 ( A1 °) 

m— 1 

In table 1 we give the values of the 's for 1 < p < 6. 
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A V) 

A\ p> 


A (V) 

AT 


1 


1 





2 





1 


3 


j 

~2 


2 


4 


3 
4 


7 

4 


5 


7 

8 


15 

8 


G 


15 
Ifi 


31 
Ifi 



where the coefficients 'are defined as 



(p) 



4"P6 



^(p) = 4 l-p _ 2 1-2 P& (P) 

= 2 3-2 P& (p) _ 2 3-2p 6 (p) + n 2 1 - 2 P&i p) 

4 P) = i2 5 - 2 P6i p) 



(A16) 



2. Two dimensions 

We start from the equation 



rn 
n 



1 p ( \ m 

rn— \ / n— 

x^r n E " )Ki3) n -^(u)' (ah) 



2=0 



By proceeding as in the case of one dimension, use of the 
algebraic relation i|Al() leads to 



where the operators Z m are defined as 

Zi = 4z/*(i) 

Z 2 = v{ii)v{i2) + v{ii)v{i 5 ) + y{ix)v{u) 
+v(i 2 )v{iz) + v(ii)v(i±) + v(h)v(ii) 
Z3 = v(ix)v(i2)v{iz) + u{ii)v{i2)v{u) 
+v{ix)v{iz)v{i ll ) + ffeM^M^) 
Z4 = v{ii)v{i2)v{h)v{u) 

(v) 

and the coefficients bm have the expressions 
4 



(A12) 



(A13) 



6l p) = 



p^ 1 / i, , 

E I = 2 p - 2 

TO =1 V "I 

J?— 1 I I II 1 .1 

E E I "" I =3(l-2P + 3P- 1 ) 

m=2 V m I n =i \ n 



P / ; ) \ " ' 

E E 



E 



2=3 \ m J n =2 \n J 1— 1 \ n 

3—1 QP I /IP— 1"> 



= 4(-l + 3 • 2P- 1 - 3 P + 4P- 1 ) 

(A14) 

By solving the system (|A12|) with respect to variables 
Z m , we can obtain the recurrence rule 



W a (i)} p =J2 A ™^ a( i)Y 



(A15) 



We note that for all p 



E^ ) = 1 



(A17) 



771=1 



In table 2 we give the values of the A$ 's for 1 < p < 8 



p 


A$ 




M 

^3 


M 

A i 


1 


1 











2 





1 








3 








1 





4 











1 


5 


3 

32 


25 
32 


35 
Ifi 


5 
2 


6 


15 

R4 


119 
R4 


75 
Ifi 


65 
16 


7 


195 
512 


1505 
512 


1799 
25fi 


175 
32 


8 


525 
1024 


3985 
1 024 


1155 
128 


1701 

256 



3. Three dimensions 



We start from the equation 

p I \ m 

m=0 \ / n=0 

>(i 2 r- n it( r V&o 

E 



X V\ 



2=0 



=0 
n—l 



rn 
n 



x > I I* 

fc=0 \ K 



(A18) 



Because of the algebraic relations <|A1[I we obtain 



[^(*)] p = 17 E b ™ Z ™ ( A19 ) 

rn — 1 
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where the operators Z m are defined as 

Z l = Qv a {i) (A20) 

Z 2 = V\V 2 + + + + v 2 v± + 

+v 2 v§ + v^vq + vav& + v*,v 6 (A21) 

Z 3 = V\ V 2 V?, + V\V 2 V{ + VlV^Vi + V 2 V^Vi + V\V 2 V^ 
+ V 2 V3 U5 + V\V A V^ + V 2 ViV^ + 
+ VlV 2 V & + V\V?,V§ + V 2 V?, V% + V\l>±l> & + V 2 l>±V & 

(A22) 

Z± = V\V 2 V?,VA + V\V 2 V?,V<i + VlV 2 ViV$ + VyV^V^V^ 
+ V\V 2 V Z V§ + VlV 2 V A V & + ViV^V^Vq 
+v 2 v?,vav$ + VlV 2 V 5 V 6 + VlV 3 V 5 V 6 + V 2 V Z V$V§ 

+v\Vav$vq + v 2 v A v$v & + v 3 v A v^v & (A23) 

Z5 = ^2^4^5 + ^1^2^3^4f6 + Vi_V 2 V z V^V & 

+vxv 2 v±v^v<o + vxv-iV^v^ + v 2 VzViV 5 v e 
Z & = v^v-sv^v^ (A24) 

and the new coefficients bm\p = 5,6) have the expres- 
sions 



p 


4(0 
Af> 


AT 


A 

^3 




^5 


A w 


1 


1 

















2 





1 














3 








1 











4 











1 








5 














1 





6 

















1 


7 


324 


/1Q 
216 


162 


245 
172 


1 7^ 
1 1 

36 


( 

2 


8 


35 
648 


1009 
1296 


2695 
648 


13811 
1296 


245 
18 


133 
18 


9 


665 
5832 


6307 
3888 


98915 
11664 


9065 
432 


10913 
432 


49 
4 


10 


245 
1296 


62167 

23328 


53375 
3888 


774575 
23328 


4165 
108 


7609 
432 



APPENDIX B: THE ENERGY MATRIX 

The energy matrix e^ d \ defined by Eq. i|26|) can be 
immediately calculated by means of the equation of mo- 
tion (|21() and the recurrence rule l12.'S[| [see Tables 1, 2, 3]. 
The matrix is defined as the matrix whose columns 
are the eigenvectors of the matrix e". In this Appendix 
we report the expressions of and for the various 
dimensions. 



m=2 \ m I n=l 




p—m Ui m—n 



5(1 - 2 P+1 + 2 • 3 p - A p + 5'P- 1 ) 



, = E P E 1 

m=2 \ m I n =l 




a. 



= -6 + 15 • 2 P + 2 3+2p - 20 • 3 p + 7 • A p - 6 • 5 P + &> 

(A25) 

By solving the system (|A19|) with respect to variables 
Z m , we can obtain the recursion rule 



c(l) 



One dimension 



-tu 2V 
-/j, 2V 
-V W - 



/ 1 2 2 1 
2 1 
yO 1 1 



(Bl) 



where the coefficients A[n are defined as 



(A26) 



Two dimensions 



^6 



^[6-36^ 
^[186^- 
^[366^- 
^[54^- 

J_t324l(p) 
6p I 5 5 
1 324 h (P) 
6p 5 °6 



2h 



(P) 



186; 

546 



(p) 
3 

(p) 



2°4 
33l(p) 
— °4 



GAP) 
5°5 



156 



(P) 



636, 



(p) 



1086^ + 1536^1 



135l(p)i 
~ u 6 J 



137 1 
10 L 



,(P)l 



1626, 



(p)i 



We note that for all p 



In table 3 we give the values of the im''s for 1 < p < 10 



(A27) 



(A28) 



,(2) 



-ft 


AV 














-n 


AV 














-V 


AV 
















AV 





8 v 


8 v 


-fv 


10V- 



( 1 4 4 2 4 

4 3 2 3 

4 2 2 2 

4 2 

\0 1 1 



(I) 4 l\ 
1 

1 

1 

i7 



(!) 3 

(i) 2 

(!) 
1 



(B2) 



(B3) 



3. Three dimensions 



1. One dimension 



,0) 





bv 


u 


u 





u 








-M 


6V 






















6V 




















-V 
























or 






















6V 





54 ^ 




203 y 
27 " 


245^ 
12 ' 




21V- 
(B4) 



(l 


6 6 


3" 


2« 


(!) 





6 s 


3 5 


2 s 


(1) 





G 4 


3 4 


2 4 


(1) 





6 3 


3 3 


2 3 


(1) 





6 2 


3 2 


2 2 


(1) 





G 


3 


2 


(1) 


Vo 


1 


1 


1 


i 



(!) i\ 
(!) 
(!) 

(!) 
(!) 
(!) 



(B5) 



/«(k) 



7 (i) 
1,1 
7 (i) 

1,2 
7 (1) 
1,3 



j(l) 
1,2 

j(l) 
1,3 

r(l) 



7 (l) 

1,3 
7 (1) 
J 2.3 

r(i) 



'2,3 ^3,3 



where 




(k) E (" 

n=l 

(m = 



L > (2)* 

2,3) 



+i) /ei + i(k) 



(m = 2,3) 



Two dimensions 



APPENDIX C: THE NORMALIZATION MATRIX 

We recall the definition of the normalization matrix 

= (^ (£i) (i^)> (d)t (j^)}) 
= l^ e i k -(R*-%)/W(k) (CI) 



/ j(2) j(2) j{2) j(2) j{2) \ 

Jii Jio Jit Ji/i J i r; 



/( 2 >(k) 



'i,i 

r(2) 

1,2 
r(2) 

1,3 
r(2) 

1,4 
r(2) 

1,5 



'1,2 
r(2) 

1,3 
r(2) 

1,4 
r(2) 

1,5 
r(2) 

2,5 



'1,3 
r(2) 

1,4 
r(2) 

1,5 
r(2) 

2,5 
r(2) 
J 3,5 



1,4 
r(2) 
1,5 
r(2) 
2,5 
r(2) 
J 3,5 
r(2) 
J 4,5 



1,5 
r(2) 
2,5 
r(2) 
J 3,5 
r(2) 
4,5 
r(2) 
J 5,5 



where 



It is straightforward to see that use of the hermiticity 
condition 1431 ) leads to the fact that we have to calculate 

only the matrix elements /{^(k) (m = 1, 2, • • • 2d + 1). 
The calculations of these is very easy when one observes 
the following anticommutating rule 

{c(i,t)[i^(i)]»,ct(j,t)} = *u[i^(i)]f 

- E (-i) n raaT*=r( P Uy<i,i)KWF-"ctat) 

n=l \ n J 

(C2) 

By taking the expectation value of l|C2|) we obtain in 
momentum space 



4 2 i(k) = l 




^(k) = «(•»" D-a(k) E (- 

n=l 






2,-"5) 






4>) = E ^ m+3) A (2 i +1 (k) 

n=l 


(m = 2 



Three dimensions 



m — 1 

lg,(k) = K^-c^k) £(-!)' 



1 



(2d)™- 1 

with the definitions 

c^ a = U^ a (i)i}j {d ^{i) 




(C3) 



«W = <K«F) (C4) 



J< 3 )(k) = 



r(3) 

1,1 
r(3) 
J l,2 
r(3) 

1,3 
r(3) 

1,4 
A3) 

1,5 
r(3) 

1,6 

\ / (3) 



r(3) 

1,2 
r(3) 

1,3 
r(3) 

1,4 
A3) 

1,5 
r(3) 

1,6 
r(3) 
i l,7 



r(3) r(3) 



1,3 
r(3) 

1,4 
r(3) 
i l,5 
r(3) 
J l,6 
r (3) 

-(3) 



1,4 
A3) 

1,5 
A3) 

1,6 
r(3) 

1,7 

-(3) 



A3) 
i l,5 
r(3) 

1,6 
r(3) 

1,7 
r(3) 
J 2,7 
r(3) 



(3) r(3) \ 

1 F. Jl 7 



1,6 
A3) 
J l,7 
r(3) 
J 2,7 
r(3) 
J 3,7 

(3) 



'2,7 -"3,7 M,7 



'5,7 



r(3) A3) A3) A3) A3) 
'2,7 J 3,7 J 4,7 J 5,7 J 6,7 



1,7 
A3) 
'2,7 
r(3) 
'3,7 
r(3) 
'4,7 
A3) 



'2,7 -"3,7 -"4,7 2 5,7 
A3) A3) A 3 ) A3) 



A3) . 
J 7,7 / 
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where 

/g(k) = l 
Ok) 



(3) n ^ _ (m-1) 



m — 1 

' E ' 

n=l 



TO — 1 

n 







/ 1 


2- 1 


2-2 


2-3 


2- 4 \ 






2- 1 


2- 2 


2-3 


2- 4 


2- 5 


1) (a) „-l 


^ = S 3 


2- 2 


2-3 




2- 5 


2- 6 


2, •••7) 




2- 3 


2 -4 


2- 5 


2-6 


2- 7 




V2- 4 


2- 5 


2- 6 


2- 7 


2- 8 / 



(D5) 



4 3) 7 (k) = E Ai m+5 )4 3 i +1 (k) 

n=l 



im 



■7) 
(CIO) 



APPENDIX D: THE SPECTRAL MATRICES 

The spectral density matrices cr^' n ^(k) can be imme- 
diately calculated by means of the knowledge of the ma- 
trices tt^ and /(<*) through Eq. J33J). 



One dimension 



a« = Si 



r(3) 





(1 





Si 



















1 




!! 


1 






1 



r(2) 



1 



-1 



V 2 



2- 1 2-2 2-3 
2 2-3 2- 4 



(Dl) 



r(5) - 



/ 1 

(!)" 


-1 


(r 

(!)" 




-2 




(!)" 
(!)" 


-3 


(!)" 
(!)" 


V(i)~ 


-4 


(!)~ 5 


/1 1 


1 


1 i\ 


1 1 


1 


1 1 


1 1 


1 


1 1 


1 1 


1 


1 1 


U 1 


1 


1 1/ 



(!) 
(!) 
(!) 



(!) 



-3 

-4 
-5 

-6 
7 



(§r:\ 

(!)" 
(!)" 
(!) 



1D6) 



(!P (!)"' (i)~ 8 / 



(D7) 



with 



51 = 7i,i - |(25Ji, a - 70/i,3 + 80/i, 4 - 32J llB ) 

5 2 = f (3/i,2 - 13/i,3 + 18/i,4 - 8/i, 5 ) 

5 3 - -4(3/!, 2 - 19/i,3 + 32/i,4 - 16/1,5) 
S 4 -f (/i, 2 - 7/i,3 + 14/i, 4 -8/ 1 , 5 ) 

S 5 = -| (3/i,2 - 22/i, 3 + 48/1,4 - 32/ li5 ) 



Three dimensions 



(D8) 



where 



51 — /i,i — 3/i,2 + 2/1,3 

5 2 =4(/!, 2 -/l, 3 ) 

53 = — /l 9 + 2/i.3 



Two dimensions 



r M = 



r (2) 



/1 





o\ 





























\o 





0/ 


/ 1 


4- 


-1 4-2 


4- 1 


4- 


-2 4 -3 


4-2 


4- 


-3 4-4 


4-3 


4- 


-4 4-5 


^4- 4 


4- 


-5 4-6 



4-3 
4-4 
4-5 
4-6 

4" 7 



4- 4 \ 
4-5 

4-6 

4- 7 
4 



-8 



(D2) 



Si 



r(2) = 



( I ••• o\ 

••• 

\o ••• 0/ 

/ 1 6" 1 ■ 

6- 1 6-2 ■ 

\6- 6 6- 7 ■ 



6- 6 \ 
6" 7 



(D3) 


= S 3 


/ 1 

3- 1 


3- 1 • 
3-2 • 


■ 3- 6 \ 

■ 3~ 7 








3- 7 • 


' 3_12 ) 






/ 1 


2" 1 • 


• 2- 6 \ 


(D4) 


= £4 


2- 1 


2-2 • 


• 2- 7 






\ 2- 6 


2- 7 • 


• 2- 12 y 



(D9) 



(D10) 



(Dll) 



(D12) 
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where 



(§)" ••• (ir!\ 



r(6) 



r(7) - 



(!) (1) ••• 
V(l)" 6 (I)" 7 ■•■ 

/ 1 x (§)■; •■ 

(!) (!) •■ 



V(ir 6 (!)■ 
(l 1 ... l\ 

1 1 ... 1 



1'" 

2) „ 

3> -' 

2 / 



-12 



(D13) 



(ir; \ 

(!)" 7 

(ir 12 y 



(D14) 



1 1 



(D15) 



here 



E - V r(p)/ (3) 

— -"m J l 



(v) 

The coefficients are given in Table 4. 



(D16) 
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^6 




1 


1 


147 

in 


406 
5 


441 
2 


315 


1134 
5 


324 
5 


2 





36 


1566 
5 


1044 


-1674 


1296 


1944 
5 


3 





-45 


1053 
2 


4149 
2 


3699 


-3078 


972 


4 





40 


-508 


2232 


-4356 


3888 


-1296 


5 





45 
2 


297 


2763 
2 


2889 


-2754 


972 


G 





36 
5 


486 
5 


468 


-1026 


5184 
5 


1944 
5 


7 





-1 


137 
10 


135 

2 


153 


-162 


324 
5 



APPENDIX E: RELATIONS BETWEEN THE 
ISING AND SPINLESS MODEL 

In this Appendix, we want to recall the main results of 
the Ising model and establish the relations between the 
two models. For an infinite chain the simplest method is 
the use of the transfer matrix method. The details of cal- 
culations are well known and can be found in many text- 
books. For example, we refer the reader to Refs42iSS^i. 
The magnetization per site is given by 



(S(i)) = 



sinh(/3/i) 



^/sinh 2 (/3/i) +e~ 4 P J 



(El) 



The two-point correlation function (S(i)S(i + j)) has the 
expression 

(S(i)S(t + j)) = <S(z)) 2 + (1 - (S(i)> V (E2) 



P 



(2) 



(E3) 



j^ 1 ' and are the eigenvalues of the transfer matrix 

v(i) - P 0J 



7 (2) = e /3J 



cosh(/3/i) + y/sinh 2 ^) +e~^J 
sh((3h) - ^smh 2 (/3h) + e~ 4 P J 



(E4) 



The three-point correlation function is given by^ 
(S(t)S(z+j)S(i+j+r)) = (S(i)f 

The relations between the Ising and fcrmionic models 
are 



h= ^(fi-Vd) 
J = -\Vd 



(E6) 



(E7) 



= j[l + 2 (S(i)) + (S(i)S(i + a))} (E8) 



n {2) = I [3 + 4 (5(i)> + (S(i)S(i + 2a))} (E9) 



A (2) = ± [3 + 7{S (i))+4(S(i)S(i + a)) 
16 

+ (S(i)S{i + 2a)) + (S(i)S(i + a)S(i + 2a))} (E10) 

By recalling l|116|) and by means of (|E7|) , the magneti- 
zation in the fermionic model has the expression 



<5(i)) = (l-T a ) 1 



v Ti - 2T X T 2 + 2T 2 2 

By oserving that T\ and T2 can be expressed as 
' £*\ _ 2e- 2 ? h 



(Ell) 



Ti = 1 - tanh I if 
T 2 = 1 - tanh 



2 



-2/3h_|_ e -4/3J 

1 - tanh(^/i) 



e 2 ^" + l 



(E12) 

it is straigtforward to see that (|E11(I is the same as (|E1J| . 
In the fermionic model, by means of l|142f> . we have 



MMi+j))-* 2 (S(i)S(i+j))-(S(i))' 



1 - (S(i)y 



= %p 



(E13) 
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where the parameter p is expressed in terms of v and 
A« by means of (|138fl . By using l|113|) and (|116|l . and 
by recalling (|115Jl and l|E12jl . it is easy to see that the 
expression of p, given by (|138l) is exactly equal to the ex- 
pression i|E3() . Then, the two-point correlation function 
of the fermionic model exactly agree with the expression 



(|E2|) of the Ising model. The parameters k^ 2 \ and 
A*- 2 ) can be calculated in the fermionic model by putting 
Ullfijl into l)112() -( ITl"4*)l . and in the Ising model by means 
of HE8(l - i|E10() . After lengthy, but straightforward, calcu- 
lations, using the relations i|115l) and (|E12|) . it is possible 
to show that there is an exact agreement. 
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